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We derive several bounds on fidelity between quantum states. In particular we show that 
fidelity is bounded from above by a simple to compute quantity we call super— fidelity. 
It is analogous to another quantity called sub-fidelity. For any two states of a two- 
dimensional quantum system (A^ = 2) all three quantities coincide. We demonstrate 
that sub- and super— fidelity are concave functions. We also show that super-fidelity is 
super-multiplicative while sub-fidelity is sub-multiplicative and design feasible schemes 
to measure these quantities in an experiment. Super-fidelity can be used to define a 
distance between quantum states. With respect to this metric the set of quantum states 
forms a part of a — 1 dimensional hypersphere. 
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1 Introduction 

By processing quantum information we wish to transform a quantum state in a controlled way. 
Taking into account inevitable interaction with an environment and possible imperfection of 
real dynamics it is then crucial to characterize quantitatively, to what extend a given quantum 
state gets close to its target. For this purpose one often uses Gdelity [I], here denoted by F. 
That quantity has also been called transition probability [5| : Operationally it is the maximal 
success probability of changing a state to another one by a measurement in a larger quantum 
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system. If both quantum states are pure, fidelity is the squared overlap between them. 

In the general case fidelity between any two mixed states is the function of the trace 
norm of the product of their square roots. Thus analytical evaluation of fidelity, or its direct 
experimental measurement becomes a cumbersome task. Hence there is a need for other 
quantities, which bound fidelity and are easier to compute and measure. 

The aim of this work is to present some bounds for fidelity and to develop experimental 
schemes to estimate it for an arbitrary pair of mixed quantum states. In particular we find an 
upper bound for fidelity by a simple quantity which is the function of purity of both states and 
the trace of their product. Since it possesses some nice algebraic properties we believe it may 
become useful in future research and propose to call it super-Hdelity. In a sense it is a quantity 
complementary to the one forming the lower bound proved in [3], and we tend to call sub- 
fidelity. For any two one-qubit states all three quantities coincide. Fidelity is well known to 
be multiplicative with respect to the tensor product. In this work we prove that super-fidelity 
is concave and super-multiplicative, while sub-fidelity is concave and sub-multiplicative. 

Fidelity can be used to define the Bures distance between quantum states and the Bures 
angle. As shown by Uhlmann in 4J the Bures geometry of the set of one-qubit states {N — 2), 
is equivalent to a three-dimensional hemisphere ^S^. The set of density operators, Qn, 
becomes the space of non-constant curvature by the Bures metric for iV > 3, [5]. 

We construct distance and angle analogous to the Bures distance out of super-fidelity in a 
similar way. With respect to this metric the set ft jv forms a fragment of a N'^ — 1 dimensional 
hypersphere with the maximally mixed state :— H/N at the pole. A linear function of 
super-fidelity was earlier used by Chen et al. (Gj to analyze the set of mixed quantum states 
and demonstrate its hyperbolic interpretation. 

This paper is organized as follows. In Section II the definition and basic properties of 
fidelity are reviewed. Sections III and IV are devoted to bounds on fidelity. In Section V 
we define sup- and super-fidelity and investigate their properties. Experimental schemes 
designed to measure these quantities are presented in section VI. In Section VII we analyze 
the geometry of the set of quantum states induced by the distance derived by super-fidelity. 
Concluding remarks are followed by appendices, in which we prove necessary lemmas and 
present the collection of useful algebraic facts. 

2 Fidelity betvifeen quantum states 

Consider an N- dimensional Hilbert space Tiiy. A linear operator p : TLn "Hn represents 
quantum state if it is Hermitian, semipositive, p = > 0, and normalized, trp ~ 1. Let Vl^ 
denote the set of all mixed quantum states of size N . 

Fidelity between quantum states pi and p2 is defined as [5] [T] , 

F(pi,p2) = (trlVpT^/^l)' = WpI'^pI'Yi. (1) 

where 1 1 • 1 1 1 is Schatten 1-norm (trace norm) , 

Pill =tr|A| :=tr\/Z4t. (2) 

Alternatively, the trace norm of an operator can be expressed as the sum of its singular values, 
||A||i = X)r=i Here (Ji{A) is equal to the square root of the corresponding eigenvalue 

of the positive matrix AA^ - see e.g. [7]. 
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There are different uses of the name fidelity. In pj the older notion transition probability 
has been renamed fidelity by Jozsa. In [8] ^/F has been called fidelity, while [5] uses Jozsa's 
notion, and to the latter convention we shall stick in calling fidelity the expression in Eq. ([1]) . 

For pure states the definition ^ is reduced to the transition probability. If one state is 
pure, pi — \ip){'ip\, then F{pi,p2) — (■(/'|p2 IV')- Hence for any two pure states their fidelity is 
equal to their squared overlap, F{il;,(t)) = KV'I*/')!^ =: k- 

Fidelity enjoys several important properties [21 UHl [HI lEl H] 7 which can also be proved on 
state spaces of unital C*-algebras. Some of them are: 

i) Bounds: < F{pi,p2) < 1. Furthermore F{pi,p2) = 1 iff pi = p2, while F{pi,p2) — 
iff supp(pi) _L supp(/32)- 

ii) Symmetry: F{pi, P2) = F{p2, pi). 

iii) Unitary invariance: F{pi,p2) — F{UpiW ,Up2U^), for any unitary operator U. 

iv) Concavity: F{p, api + (1 - a)p2) > aF{p, pi) + (1 - a)F(p, P2), for a e [0, 1]. 

v) Multiplicativity: F{pi ® p2,P3® Pi) = F{pi, p3)F{p2, pi)- 



vi) Joint concavity: ^/F{apl + {l-a)p2, ap[ + {l-a)p2) > avF{pi, p[) + {l-a)VF{p2, p'2), 
for a e [0, 1]. 

For further analysis of fidelity properties it is instructive to work with eigenvalues of a ma- 
/1/21/2 

trix y Pi P2P1 ■ Let us denote them by A^, i — 1, . . . , N . This matrix is positive so its 
eigenvalues and singular values coincide. Unless otherwise stated, we tacitly assume that 
Ai > A2 > • • • > Aat. The root fidelity reads 

N 

\/F{pi,p2) ^ tr^J ^/p{p2^/pi ^^K- (3) 

1=1 

Squaring this equation one obtains a compact expression for fidelity, 

2 



F{pi,p2)=\Y^^i\ =trpip2 + 2^A,Aj, 

\i=l / i<j 



(4) 



where we have taken into account that trpip2 — '^'^\fp\P2^fp\ = St=i'^f- "^^^ matrix 
y/piP2^/pi is similar to pip2 and they share the same set of eigenvalues. 



3 Bounds for fidelity 

We shall need some further algebraic definitions. For any matrix X of size N with a set of 
eigenvalues {Ai, . . . , Ajv} we define elementary symmetric functions Sm{X) as the elementary 
symmetric function of its eigenvalues [T71 Def. 1.2.9]. For instance, the second and third 
elementary symmetric functions read 

S2{X) = ^A.A„ (5) 

i<j 

i<j<k 
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For any matrix of rank r the highest non- vanishing symmetric function reads Sr{X) = Y[l=i ^i- 
In the generic case r = N we have sn{X) = det(X). 

In this section we shall list several bounds for fidelity, some of which are well known in 
the literature. Let us start by stating a simple result, 

-F(pi,P2) < trpitrp2, (7) 

which follows directly from Fact[l] (see Appendix A) if we set = 1/2. This fact implies the 
property F{pi,p2) <1. 

Expression (|4]) implies the following lower bound 

trpiP2 < F{pi,p2) < A^tr|pip2|. (8) 
To get the upper bound we use Fact[2](see Appendix A) and set — 1/2 to obtain | l-ypiy^l |i < 

N\\PlP2\\l- 

Let us now denote the spectra of the states pi and p2, by vectors p and q, respectively. 
The fidelity between them is then bounded by the classical fidelity between diagonal density 
matrices [13] 

F{p\q^)<Fip,,p2)<F{p\q^), (9) 

where the arrows up (down) indicate that the eigenvalues are put in the nondecreasing (non- 
increasing) order. 

The lower bound in ([8]) can be improved, since the following result is true [3] 

F{pi,P2) > trpip2 + V2^/ (trpiP2)^ - tlpiP2PlP2- (10) 

The above inequality is saturated for any pair of one-qubit states. Furthermore, the above 
inequality is an equality if the rank of pip2 does not exceed two. On the other hand, the 
inequality is strict if that rank is larger than two — see Appendix E. For completeness we 
present the simple proof of inequality (jlOp in Appendix B. 

Another lower bound is obtained if the rank of pip2 is exactly r. If Sr denotes the r*** 
elementary symmetric function then 

F{pi,P2) > tr/9i/92 + r{r - I) {/ Sr{piP2)- (11) 

This bound is proved in Appendix C. If both states are generic, i.e. if they are of the maximal 
rank the above formula reads 

F{pi,p2) >l^piP2 + N{N ~1) Vdetpi detp2. (12) 

The key result of this paper consist in the following upper bound, in a sense complementary 

to uni). 

Theorem 1 For any density matrices pi and p2 we have 

F{pi,P2)< trpip2 + V'(l-trp2)(l_trp2). (13) 
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Before presenting the proof in the subsequent section let us first note that the bound is 
saturated if at least one of the states is pure. Furthermore, an equality holds for any two 
mixed states of size N = 2. To show this property observe that in this case the sum in 
dH) consists of a single term 2A1A2 = 2-\/det(pip2) — ^2det(pi) Y/2det(/92)- Since for any 
one-qubit state one has 2det(p) = 1 — trp^ an equality in (fTH]) follows. This fact was already 
known to Hiibner [14j . In a similar way we treat the more general case of = 3 in Appendix 
IeI for which some other equations for fidelity are derived. 

4 Proof of the main upper bound 

The notion of the second symmetric function ([5|) allows us to write the expression 

[itTX)^~{trX')]^2s2{X). (14) 

Note that if X has nonnegative eigenvalues then S2{X) > 0. 
Using (fT4|) we can rewrite fidelity 



Fipi,p2) = trpip2 + 2s2 P^^ P2py'^ j , (15) 
and 

^(1 - trp2)(l _ trpj) - 2^S2{pi)s2ip2). (16) 
Thus the Theorem [T] can be equally expressed as an inequality 

< V^2{PI)S2{P2)- (17) 

The proof of (|17p is decomposed into two Lemmas, the proof of which can be found in 
Appendix |DJ 

Lemma 1 For given density matrices pi, p2 with eigenvalues pi, . . . ,p„ and 51, . . . , (?„ respec- 
tively 

S2 (y\J p\''^ P2p\''^^ < S2 (Vdiag(p)diag(q)) , (18) 

where diag(p) and diag(q) denote diagonal matrices with entries on diagonal pi, . . . ,p„ and 
<?i , . . . , (7,1 respectively. 

Lemma 2 With notation as in LemmaUl we have 



S2 



S2 



(^v/diag(p)diag((7)) < Vs2(diag(p))s2(diag(g)) = \/s2(pi) 52(^2)- (19) 



Proof of Theorem [TJ For given density matrices pi and p2 with eigenvalues pi, . . . ,P' 
and gi , . . . , g„ respectively. We denote diagonal matrices with entries on diagonal pi , . . . , p„ 
and qi, . . . ,qn as diag(p) and diag((7) respectively. 



F{pi,P2) = trpip2 + 2S2 (^^J p\''^ P2p\''^^ 



< trpip2 + 2s2 ( V diag(p)diag(q) 



< trpip2 + 2v/s2(pi)s2(/32) = trpip2 + Y (1 - trp^)(l - trp^). 

Making use of Lemma [1] for the first inequality and of Lemma [5] for the second one we arrive 
at the inequality (|13p. □. 
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5 Sub and super fidelity and their properties 
5. 1 Definition and basic facts 

We shall start this section with a general definition. For any two hermitian operators A and 
B let us define two quantities 

E{A, B) = tiAB + ^2[{iT: ABY -tiABAB], (20) 
G{A,B) = tvAB + V(trA)2 _ trAVltr^)^ " tr-B2. (21) 

For any two density operators their traces are equal to unity, so E{pi,p2) and G{pi, P2) 
have lower bound pUj) and upper bound p^ . respectively. Thus both universal bounds for 
fidelity can be rewritten as 

E{pi,p2)<F{p^,P2)<G{pi,p2). (22) 

Note that both bounds require the evaluation of three traces only, so they are easier to 
compute than the original fidelity As shown in Section [3] for iV = 2 all three quantities are 
equal, so we propose to call E{pi,p2) and G{pi,p2) as sub- and super-Gdelity. These names 
are additionally motivated by the following appealing properties: 

i') Bounds: < E{pi,p2) < 1 and < G{pi,p2) < 1. 
ii') Symmetry: E{pi,p2) = E{p2,pi) and G{pi,p2) = G{p2,pi)- 

iii') Unitary invariance: E{pi,p2) = E{U piU\U p2U^) andG(pi,p2) = G{U piU\U p2U'^), 
for any unitary operator U . 

iv') Concavity: 

Proposition 1 Sub- and super-fidelity are concave, that is for A, B,C, D € Qj^ and 
a S [0, 1] we have 

EiA,aB + {l-a)C) > aE{A, B) + {1 - a)E{A, C), (23) 
G{A,aB + {l-a)C) > aG{A, B) + {1 - a)G{A,C). (24) 

v') Properties of the tensor product: 

Proposition 2 Super-fidelity is super-multiplicative, that is for A, B,C, D £ 

G{A(E)B,C(g)D)>G{A,C)G{B,D), (25) 

while 



Proposition 3 Sub-fidelity is sub -multiplicative, that is for A,B,C,D € f2jv 

E{A® B,C ® D) < E{A,C)E{B,D). (26) 
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Properties i'), ii') and iii') follow from the properties of trAB and definitions ([20|l and ([2T|) . 
In this section we prove properties iv') and v'). 

Proof of Proposition [H The definitions ^ and ^ can be rewritten in terms of the 
aforementioned elementary symmetric functions ([5|) using relation (|14p . 

E{A,B) = tYAB + 2^/s2{AB), (27) 
G{A,B) = tYAB + 2^s2{A)s2{B). (28) 

In general r^^ root of Sr is concave on the cone of positive operators [41]. This implies 
concavity of G directly. To get concavity of E we can replace matrix AB by the similar 
matrix A^^^BA^^^ which is positive. Using the concavity of \/ S2{A^/'^BA^/'^) we obtain the 
result. □. 

Proof of Proposition [2j First we note that super-fidelity is not multiplicative. As an 
example we can take 

--{I D^c^c i).o.{i °), ,2.) 

in which case we have 

^ = G{A ®B,C®D)> G(A, G)G{B, D) = 0. (30) 
To prove the proposition we write 

G{A ®B,G®D)^ iiACiiBD + ^(1 - trA2trB2)(l - IiCHtD'^), (31) 

and 

G(A, C)G{B, D) = (trAC + ^{l ~ trA'^){l - trC^)^ (trBD + {I - tT:B'^){l - trD'^)^ . 

(32) 

Denoting tiA^ = a, tiB^ = f3, trC^ — 7 and tr_D^ = 5 we have to show that 

\/(l-a/3)(l-7<5) > ti-ACy/(l - /3)(1 - <5) + trBDy/{l - a){l - 7) 

+ v/(l-«)(l-7)(l-/?)(l-'^)- 
Now from Fact [5] with a — 2 (see Appendix B) one has 

tr AC < \/tryl2trG2 = (33) 

and 

ti-BD < VtrBHiD^ = y^. (34) 

Thus it is enough to show that 

V(l-«/3)(l-7'^) > V ^Va -p)il~6) + V^V(1 - a)(l - 7) 

+ ^{l-a)il-j)il-m-S). (35) 

We define two vectors 

X = V^\/r^^ and y = VSVT^j_ . (36) 
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Note that 

{X\Y) = - m - <5) + - - 7) + V(l - - 7)(1 - - S) (37) 

and 

{X\X) = (1 - aP) and {Y\Y) = (1 - 7(5). (38) 
Now by combining (j38p with (|37p and using Cauchy-Schwarz inequaUty 

V{X\X){Y\Y)>(X\Y), (39) 

we obtain ([55]) . □. 

Proof of Proposition [3l To show sub-multiplicativity of sub-fidelity we write the defini- 
tion (PO)) for a tensor product, 

E{A^B,C(giD) = tr[{A(^ B){C ® D)] 

+ y^2[(tr(yl (g) B){C D)Y - tr{A ® B)(C ^D){A(» B){C ® D)] 
= iiACiiBD + ^2[{tiACtTBDY - trACACtvBDBD]. 

The product of two sub-fidelities reads 

E{A, C)E{B, D) = (trAC + y/2[{tT:ACY - tr ACAC] ) (trSZ) + v/2[(trBi:))2 - trBLiBZ)]) 
= trACtrBZ) + ivAC ^/2[{iT:BD f - liBDBD] + trSD V2[(tr^C')^ - ItACAC] 
+ y/2[{ir ACf - trACAC]\/2[(trBD)2 - tiBDBD]. 

For short wc denote 

a = tvAC, a = tr AC AC, 
P = trSD, b = tvBDBD. 

We have > a and /S^ > 6. By rewriting above expressions in the new notation we 
obtain 

E{A®B,C®D) = a/3 + ^/2[a'^(3^ - ah] 

and 

E{A,C)E{B,D) ^ a(3 + a^2[l3^ -h]+ ii^j2{a^ - a] + ^2\o? - a\^j2\iS^ - h\. 
Now we write 

= v/2[(a2 - a)(/32 - 5) + 6(0^2 - a) + a(/32 - h)\. 
Making use of subadditivity of square root we obtain 

V2[a2^2 _ < ^2(a2 „ a)(/32 - 6) + ^26(^2 - a) + V2a(/32 - 6)]. 
Because 2 < A, a < a'^ and 6 < /32 we get 

V2[a2/32 - a5] < ^4(^2 _ a)(-^2 _ 5) _^ /3^2(a2 _ a) + a^2(/32 - b)]. 
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And as a result we obtain the desired inequality 

E{A(g, B,C(S D) < E{A,C)E{B,D). 

□ . 

For any pair of Hermitian operators Xi and X2 let us now define a quadratic Lorentz form 

(Xi,X2)l [(trXi)(trX2) - trXi^a]. (40) 

To find out the motivation standing behind this name let us expand a Hermitian operator 
X in an operator basis, X — J2f=o^ '^j^j- We assume that the basis is orthogonal, tvHjHk — 
Sjk, the first operator is proportional to identity, Hq = 1/ x/iV, and all other operators Hj are 
traceless. Then the form (|40p gives 

(X, X), ^ {irX f trX' = ^^^a^ - (^l) 

which is of Minkowski-Lorentz type. By the help of this notion, super-fidelity can be written 
as 

G{A, B) ^ irAB + [(A, A)l{B, B)lY'\ (42) 

while sub-fidelity reads 

E{A,B) =iiAB + [2{AB,AB)lY/'^. (43) 

The forward cone with respect to the form (|40p is given by operators X satisfying 

{X, X)l > and trX > 0. (44) 

Since the density matrices are normalized, trp = 1, the form {p,p)l is non-negative. 
For a Lorentz form any two forward directed Hermitian matrices A and B satisfy 

{A,A)l{B,B)l<[{A,B)l?. (45) 

Substituting this bound into expression (|42p we arrive at an upper bound for super-fidelity 

G{A, B) < iiAB + (trA)(trB) - {tiAB) = (tryl)(trB). (46) 

For the case of normalized density matrices, ivp = 1 we get G{pi,p2) < 1. 
Using the bound (jlip for density operators we introduce a third quantity 

E'{A, B) = tvAB + r{r - 1) ^Sr{AB), (47) 

where r is the rank of matrix AB. Note that for r = 2 this formula is reduced to an expression 
(P7)) for sub-fidelity, hence in this case E' ~ E. 

Since {sr{X)Y/^ is concave for density operators we infer that the quantity E'{A,B), 
defined in equation (|47p . is separately concave in A and in B. 
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5.2 Examples and classical analogues 

To observe sub- and super-fidelity in action consider a family of mixed states 

Pa^a\ij){^\ + {l-a)llN, 



(48) 



which interpolates between arbitrary pure state and the maximally mixed state. It is 
straightforward to compute the fidelity between the state pa and the maximally mixed state 

F{Pa, P*) = ^ [^[N -l)a+l + {N- l)VT^) ' , (49) 



as well as other bounds 

E{pa,P*) = 

E'{pa,P*) = 

G{pa,P*) = 



N 



N 



N 



V«^-l)a + l)(l-a) 



N-l 



N 
1 

iV 



These results are plotted in Fig. [T]for TV = 2, 3, 4, 5. 




(50) 
(51) 
(52) 




b) A' = 3 




Fig. 1. The comparison of sub-fidelity E, bound E' , fidelity F (solid line) and super-fidelity 
G. Each plot shows these quantities calculated for the maximally mixed state and a state l|48|l 
depending on the parameter a. For a one-qubit system, case a.) N = 2, one has E = E' = E = G. 
Note the difference between these quantities shown for A'^ = 3, 4, 5. In this case E > E' for a close 
to unity. 



For iV = 2 all these quantities coincide, and the quality of the approximation goes down 
with the system size N, as expected. Looking at the graph one could imagine that relation 
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E < E' is fulfilled. However, such an equality does not hold as we found a counter example: 
the pair of states analyzed in the figure with parameter a very close to unity. 

One may work out several other examples, for which sub- and super-fidelity are easy to 
find. Explicit formulas are simple in the case of two commuting density matrices pp and pq 
with spectra given by vectors p and q, respectively. In such a classical case these quantities 
read 



E{Pp,Pq) 

F{Pp,Pq) 
G{Pp,Pq) 





N 



N 



(53) 

(54) 
(55) 



i=l 



F 



5.3 The difference G 

In view of the inequality (|22p it is natural to ask how big the difference G — F might be. 
Since both quantities coincide if one of the states is pure, let us analyze the case of two mixed 
states living in orthogonal subspaces. 

More precisely, let us fix an even dimensionality of the Hilbert space N — 2M, and define 
two diagonal states, each supported in M dimensional space, pi = ■^diag(l, . . . , 1, 0, . . . , 0) 
and p2 = ;^diag(0, . . . , 0, 1, . . . , 1). Since they are supported by orthogonal subspaces their 
fidelity vanishes, F(pi, P2) = 0. On the other hand the definition (pij) gives their super-fidelity 



G(pi,P2) 



N 



N 



(56) 



equal in this case to the difference G— F. As expected for = 2 we get G = F = Q. However, 
for N large enough the difference G — F may become arbitrarily close to unity. 

Thus working with super-fidelity G in place of fidelity F one needs to remember that this 
approximation works fine for small systems or where at least one of the states is pure enough. 

6 On measurement methods 

6.1 Associated physical observables 

Here we shall shortly discuss possibilities of measurement of both sub- and super-fidelities 
in physical experiments. The approach below follows the techniques used in state spectrum 
estimation [20 and nonlinear entanglement detection and/or estimation which has been devel- 
oped significantly last years (see [HI |^ and references therein) . Those approaches exploited 
the properties of SWAP operator and other permutation unitary operations to get the prop- 
erties of single state rather than the relation of different states. There were little exceptions: 
one was a quantum network measurement of an overlap of the two states [20'. Here we shall 
follow the latter idea since we want to estimate the distance of two different quantum states. 
In particular we shall see that it is possible to measure these quantities with help of not 
more than two collective observables. This fact may be helpful in experimental comparison of 
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two different stationary sources of quantum states. Quite remarkably, as we sliall see below, 
with help of similar techniques, super-fidelity can be represented by only three experimental 
probabilities which makes it very friendly from an experimental point of view. 

We start by providing a simple example. First one can see that to calculate sub- and 
super-fidelity it is necessary to calculate the values of the terms of the form ti AB. Let 
A.B e Af2(C)- In this case 



A = 



ail ai2 
0.21 022 



,B = 



hi bi2 
621 622 



(57) 



and 
tvAB = tr 



011611+012621 011612 + 012622 
021611 + 022621 021612 + 022622 



011611 + 021612 + 012621 + 022622- (58) 



On the other hand this value can be calculated using SWAP gate as 



tr[SWAP{A(g)B)] = tr 



= tr 



/ 1 




V 

/ 011611 
021611 
011621 

V 021621 



00^ 
1 

1 / 

011612 
021612 
011622 
021622 



f 011611 
011621 
021611 
V 021621 

012611 
022611 
012621 
022621 



011612 
011622 
021612 
021622 

012612 ^ 
022612 
012622 
022622 J 



012611 012612 \ 

012621 012622 

022611 022612 

022621 022622 / 

= trAB. 



(59) 



(60) 



To address the question of measurability of the quantities (|20p . (PTjl let us first recall the 
corresponding permutation operators which we shall need subsequently. The first one will be 
just SWAP operator (example of which is the SWAP gate presented above) ^12 : 'Hn®'Hn 
TLn ^'Hn which is defined by the action 



Vl2\(t>l) ® IV'2) = IV'2) <^ 101 



(61) 



This is a Hermitian operator and as a such it represents an observable. It has a simple 
eigendecomposition in the form 



V12 



12 



(-) 



12 ! 



where projections Pi2 ^ onto symmetric and antisymmetric subspaces of Hn (8) TYat are 

Pi±2 = i(Il2±1^12). 



(62) 



(63) 



Below we shall omit the indices and use the notation V and P^^\ if it does not lead to 
confusion. An important property usually exploited in case of entanglement detection is that 
the formula ((60)) holds for the SWAP operator V of any dimension [23]. Apart form that 

-ij(84 -1/(84 

'I- AT ^ TV ' 

(64) 



operation we will also need a family of unitary permutation matrices V-^234 ■ ' '-n 

V123M 1^-2) <» IV-s) <^ 1^-4) = l^^(l)) ® 1^-^(2)) ® IV'^O)) ® 1^-^(4)), 



where tt represents any chosen permutation of the indices (1,2,3,4). For simplicity we shall 
drop the indices using the notation . 
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Let US define the set S of all eight permutations that do not map the sequence (1, 2, 3, 4) 
into a one having odd or even elements one after another. For instance, the permutations 
defined by the ranges (2341) or (3214) belong to S, while (2314) or (1423) do not. For a fixed 
set S' C S and some ttq € 5 we define the following observable: 



A special case is the observable W^'^°^''^° = {V'^° + V^o )/2 with ttq being just some cyclic 
permutation (cf. [21] and references therein). The choice of the permutation ttq and/or the 
subset 5' may be motivated by a specific physical situation. 

In the case of single qubit sources {N = 2) the observables ((65)) have highly degener- 
ated spectra and the corresponding eigenvectors have very symmetric forms. In particular 
the observable W^'^°^''^° has spectrum {1,-1,0} which means that its mean value requires 
probabilities of only two outcomes of incomplete von Neumann measurement. The observable 
has the spectral decomposition W^'^°^''^° = Q^^-' — Q*^^-* where support of the projector Q^^-* 
is spanned by eigenvectors {|0i) = |0000), |(/i2) = \h) = (lOlU) + |1011) + |1101) + 

|1110))/2, |(/.4) = (10011) + 10110) + llOOl) + |1100))/2, 105) = (10101) + |1010))/V2, \q^e) = 
af'^lt/ji)} while the support of the second projector Q^^^ (orthogonal to Q^^-*) corresponds to 

To illustrate how to measure the quantities E and G suppose now we can perform collective 
measurements on two and four copies of both quantum states. We plan measurements that 
allow two or four copies of analyzed states to interact. Then (cf. [20l[2Tl[22] and references 
therein) the sub- and super-fidelities can be represented in terms of averages of following 
observables, 



E{pi,p2) = trFpi «) p2 + V2[(trV^pi (g)p2)2 - trl^-S.Trop^ (g, P2 (g) «) Pa], (66) 



There are two simple but important observations to be made. The sub-fidelity E can be 
measured with help of two setups: (i) the one measuring the observable V and (ii) the second 
one measuring observable W^''^° . Each setup requires one source: setup (i) needs the source 
that creates, say, pairs pi® P2, while setup (ii) requires a source producing quadruples of the 
form, say, pi ® P2 ® Pi ® P2- 

Our scheme will work also for a worse source that produces one of the pairs (quadruples) 
{pi ® P2, P2® Pi} {{pi ® P2® Pi® P2, P2®Pi® P2®> Pi}) at random according to an unknown 
biased probability distribution, which will not affect the results of the corresponding estimate 
for sub-fidelity. 

The second observation is that the super-fidelity G can be measured with help of single 
setup, namely the one that measures observable V, but requires its application to three types 
of sources i.e. the ones creating pairs pi ®) pi, P2® P2, and, say, pi ® p2- Again, the last 
source may produce at random one of the pairs {pi ® P2i P2® Pi} and this will not affect the 
estimate for super-fidelity. 

It is very interesting to study the form of super-fidelity in terms of directly measurable 
quantities, i.e. probabilities^ since it has a simple optical implementation. Let us introduce 




(65) 




G{P1,P2) 



irVpi I® P2 + trVpi ^ pi\^l- trVp2 P2 ■ 



(67) 
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the probabilities of the projection onto the antisymmetric subspace of Ti^v <8) Hat : 

p|7^ =trP(-V,0pj, z,j = l,2. 



(68) 




(69) 



which is crucial for further discussion. Note that the super-fidelity can be represented in 
terms of only three probabilities that can be measured in a single set-up. One can perform 
a simple consistency test by checking, whether the combination of experimental probabilities 

satisfy (up to error bars) the condition ^ ~ \/ Pii ^22 ^ — 0-^ ~ otherwise one had an 
unphysical result, since super-fidelity can not be negative. Note that the probability p^-^^ ^ 
has been already measured experimentally for two copies of composite systems in context of 
entanglement detection [25 ] or estimation |26j under some assumptions about the nature of 
the sources. In subsection below we shall refer to the scheme analogous to the one utilized in 



It is interesting to note that if the state g is of d-dimensional type, then reproduction 
of sub- and super-fidelity via quantum tomography requires 2(P — 2 independent quantities 
to be estimated since each of the two states is described by — 1 real parameters. On the 
other hand, to find the quantities E and G in the way described above one requires only two 
or three independent real quantities (probabilities) to be estimated independly on how large 
the dimension d is. The price to be payed is, of course, that one must perform collective 
experiments. Preparation of reliabe setups of such experiments might be a good test for 
quantum engeneering. 

6.2 Measuring super fidelity of states representing photons polarizations 

Consider now physical setup that would compare two states of polarization of single photon 
in terms of super-fidelity G. In this case the density matrix is defined on Hilbert space 
isomorphic to where the horizontal (vertical) polarization, usually denoted by \H) {\V)) 
corresponds to the standard basis element |0) (|1)). Suppose one has memoryless sources of 
two types Si [i — 1, 2) sending photons in polarization states pi, p2 respectively. 

The experimental setup is elementary. We have sources Si, Sj, where we put either 
i = j = 1,2 (sources of the same type) or, say i = 1, j = 2 (different sources) then we have 
a beamsplitter (in equal distance to the source) and two detectors behind it (see Fig. If 
two photons form sources Si, Sj meet on the beamsplitter and the two detectors click, we 
have so-called anticoalescence event, which happens with probability [25 . Otherwise 
we deal with a coalescence result which occurs with probability = 1 — p[j . Putting 
all three probabilities of anti-coalescence into formula (|69p we reproduce the expression for 
super-fidelity. 

This seems to be the most easy experiment with two sources to perform. Such an experi- 
ment can be realized for two sources of photons engineered with help of controlled decoherence 
(in a way similar to Ref. [27j ) corresponding to two different mixed states of a qubit. 

The above scheme immediately extends to the case of states pi , p2 are defined on iV = 2"- 
dimensional Hilbert space representing polarization degrees of freedom of n photons. In this 



Ref. m 
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Di D'i 




Fig. 2. Elementary scheme with one beam-sphtter BS. The sources Si, Sj are chosen to be, in turn, 
ot the same {i = j = 1,2) and different (i = 1, j = 2) type. Single click in either of the detectors 
D, D' corresponds to projection into symmetric two-qubit subspace of photon polarization, while 
two clicks represent projection onto ono-dimensional antisymmetric (singlet) subspace. 



case the total Hilbert space is TCn — (C^)®" and the scheme reads as in Fig. [3] (compare 
[251 [^). If the probabihty with = —1 (sfc — +1) corresponds to anticoalescence 

(coalescence) on fc-th beamsplitter, i.e. it represents the probability of two clicks (one click) 
in the pair of detectors 1?^, D'f^, then the total probabilities: 

P^'= E pf^'Mf^'^..pg"^'", M = l,2 (70) 

si,S2...,s„: siS2...s„ = — 1 

are these we put into ([M)) . In the formula above we count all the cases when an odd number of 
anti-coalescence events occurs, provided that there is no photon losses during the experiment. 

D, D'i D2 D'2 Dn D'n 




Fig. 3. The scheme for measurement of super— fidelity of the states of n-photon polarizations. 
According to jTOj only the events with double clicks in odd number of detector pairs contribute 
to each of the three probabilities in the formula II69I I. 

For n = 2 this type of experiment has already been performed with two two-photon 
sources producing entangled states [25j . However, the sources were considered to provide 
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the same state on average rather than two different ones. A similar reasoning was used in 
another recent experiment, in which photon polarization and momentum degrees of freedom 
were used to estimate the concurrence where additional strong assumption about purity of 
each copies were also used |26j . In general, measurements schemes of quantities like purity, 
concurrence, sub- and super-fidelity in the collective framework like the one presented here 
requires the assumption that the sources producing states are stationary and memoryless. 
Quite remarkably this is the same assumption one makes in quantum tomography. As dis- 
cussed in [29, ,30J there may be difficulties with satisfying it in real experimental scenarios for 
instance due to classical correlations between the consecutive copies of the system. In other 
words the condition of having the global state in form (which mathematically corresponds 
to quantum de Finnetti condition [21]) may not be obeyed. This point requires more analysis, 
and leads, in general to nontrivial issues. It seems however that stationarity and memoryless 
character of the source can usually be satisfied approximately. Then the present measure- 
ment would serve (similarily like quantum tomography does, though may be in a way more 
sensitive to source correlations) as an approximative, coarse-grained-like characteristic of the 
states sources under reasonable physical assumptions. 

6.3 Quantum networks 

There is yet another method of detection of quantities that may be considered here. This is a 
method based on quantum networks. It is known that a unitary operation U acting on state 
a but controlled by a qubit in the superposition state |-|-) = jdO) + |1)) reproduces the value 
i?e (trJ7 cr) directly as a mean value of the Pauli matrix {a^) measured on the controlled qubit 
[24l [20] . This fact allows us to measure certain nonlinear functions of the state. To get trp'' 
one takes k copies of the state, a = p®^ and takes for U the operator of cyclic permutation, 
in full analogy to used in previous subsection. To measure the overlap of two matrices, 
pi, p2, one takes cr = pi ® p2 and uses the SWAP operator, U — V. 

The corresponding network is already provided explicitly in Rcf. [20 , so we shall not 
write it down here. Such a network allows one to measure all three quantities needed to 
reproduce the super-fidelity G. Indeed, the network produces directly (as mean values of CTz 
on controlled qubit) all three mean values: tvpiPj, i = 1,2 provided that states forming the 
input of the controlled part of the network are pi (gi Pj- Some alternative constructions of 
programmable networks designed to measure super-fidelity are also possible. 

Similarly, sub-fidelity E can also be estimated with a network-based experimental scheme. 
Following reasoning from Ref. |21j one constructs the following programmable quantum 
network - see Fig. |4] Depending on the program state |5'i2) as a mean value of CTz of the 
measured controlling qubit one gets 

(i) trpip2 if 1*12) = |0)|0), 

(ii) trpip2PiP2 if 1*12) = |1)|0), 

(iii) i (trpip2Pi/02 - (trpiP2)^) if 1*12) = (|0)|1) + |1)|0))/V2 i.e. if it is in BeU state. 

The last quantity up to the factor (— -j) is just the quantity that occurs under the square 
root in the formula for E. In general to estimate the sub-fidelity E we may "run" the 
first "program" and either the second or the third one. Alternatively, we may run all three 
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|«'12) 

|0> 
|0) 

Pi 

P2 
Pi 
P2 



X 



X 



^ e 



-{a\ [W 



Fig. 4. Example of programmable network allowing to measure in particular the quantities 
tr(pip2PiP2) and ^tr(pip2PiP2) — tr(pip2)^. The state |\I'i2) represents the program (see the 
main text) . Symbol M corresponds to the measurement of Pauli matrix ■ The last phase gate 
usually is chosen to be a = unless the interferometric picture with visibility is needed (cf. Ref. 
[20]). See [8] for the description of quantum gates used in this circuit. 



programs and use the data to verify the accuracy of the experiment by comparing the two 
partially independent estimates of E obtained in that way. It is easy tu see, that the same 
network can be used to estimate super-fidelity G if one puts as an input pi® Pi® Pj (8" Pj, 
J = 1,2. 

7 Distances and geometry of the space of states 
7.1 Hubert- Schmidt distance and flat geometry 

The geometry of the space of quantum states depends on the metric used [32l [33l [SH [9] . The 
set rijv of mixed states of size N reveals the Euclidean (flat) geometry if it is analyzed with 
respect to the Hilbert-Schmidt distance, 

Dhs{pi,P2) - ^ti[{pi- P2?]. (71) 



To demonstrate this property let us first concentrate on the simplest case, N ~ 2. Making 
use of the notion of a coherence vector f any state of a qubit can be written in the Bloch 
representation 

p=l+f-X. (72) 

Here A denotes the vector of three rescaled traceless Pauli matrices {ct^, ct^, ct^}/ V2, which are 
orthogonal in the sense of the Hilbert-Schmidt scalar product, (A'^jA™) = tr(A'^)tA™ = 6km.- 
Together with A° = I/\/2 they form an orthonormal basis in the space of complex density 
matrices of size two. Due to Hermiticity of p the three-dimensional Bloch vector r is real. 
Positivity condition implies |r| < 1/V2 — R2 with equality for pure states, which form 
the Bloch sphere of radius R2. Representation ([7^ implies that for any state of a qubit 
trp2 = l/2+|T|2. 

Consider two arbitrary density matrices and express their difference pi — p2 in the Bloch 
form. The entries of this difference consist of the differences between components of both 
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Bloch vectors fi and f^. Therefore Hilbert-Schmidt distance induces the flat (Euclidean) 
geometry of Q2, 

(Pfi , Pf2 ) = -Db (n , 7^2 ) , (73) 

where -D_b is the Euclidean distance between both Bloch vectors in ]R^. 

It is worth to add that expression ((73|l holds for an arbitrary TV. In this case r is a real 
vector with iV^ — 1 components, while the vector A = {^''}k=i^ ((72|) denotes the set of 
N'^ — 1 traceless generators of the group SU{N). Positivity of p implies that the length of 
the Bloch vector is limited by 



N -1 

For N = 2 the condition |r| < R2 is sufficient to imply that the corresponding matrix is 
positive and represents a state, while for TV > 3 it is only a necessary condition [9J. This is 
related to the fact that with respect to the flat, H-S geometry the set ^2 forms a full 3-ball, 
while for larger N the set ^In forms a convex subset of the {N"^ — l)-dimensional ball of radius 
Rn centered at = 1/N. 

7.2 Bures distance and the geometry it induces 

The notion of fidelity, introduced in ([1]), can be used to define the Bures distance [35l [2] 



DFipi,P2) = ^2-2Vi^(pi,^. (75) 
or the Bures length [35] (later called angle in [5]), 



1 



D'p{pi,p2) arccosv/F(pi,p2) = -arccos(^2i^(pi, P2) - 1 j ■ (76) 

For any pair of pure states the Bures length coincides with their Fubini-Study distance, 
D'p{pi,,p<i>) = dFs{\ip), I'/')) = arccos|(V'|0)|. 

The Bures metric is distinguished by its rather special properties: it is a Riemannian, 
monotone metric f37^, Fisher adjusted metric [33], closely related to the statistical distance 
[32]. 

It is not difficult to describe the geometry of the set of mixed states of a single qubit 
induced by the Bures metric. Consider a mixed state p (E ^2 and its transformation proposed 
mi 

P{x, y, z) ^ (^x, y, z, t = y i - a;2 - 2/2 _ ^2'^ ^ (77) 

It blows up the Bloch ball B'^ of radius R2 = into a hyper- hemisphere ^S^ of the same 

radius. The original variables (cc, y, z) denote the parameters of the state in the Bloch vector 
representation. The auxiliary variable reads t = -^1/2 — |tP in terms of the Bloch vector, so 
that t^ + trp^ = 1. The maximally mixed state p, — (0,0,0), is mapped into a hyper-pole. 
It is equally distant from all pure states located at the hyper-equator S*^, which form the 
boundary of ^2- 

Any state p is uniquely represented by an 'extended Bloch vector', v = {x, y, z, t) of length 
i?2- The auxiliary variable reads t = ■\/l/2 — |rp in terms of the Bloch vector, so that 
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t'^ + trp^ — 1. Consider two states pi and p2, described by two vectors vi and V2 G IR-^, 
which form the angle z9. Since for any one-qubit states the bound (|13p becomes an equahty, 
we see that fidelity between them reads F(pi,p2) = trpip2 + — ^/'^ + Ti ■ T2 + ^1^2- 

This can be rewritten with the use of extended vectors Vi and the angle between them, 
F ^1/2 + vi ■ V2 ^ 1/2 + cos-d. Since i?i = 1/2 we find 

■d = aiccos{2F^l) =2D'p{pi,p2), (78) 

which shows that the Bures length ((76)) between any two mixed states is proportional to the 
Riemannian distance between the corresponding points at the Uhlmann hemisphere. 

Making use of the fidelity F one can also define other distances in the space of quantum 
states. For instance, Gilchrist et al. [38 have shown that the root infidelity 

C{pi , P2) = V^-Fipi,P2) (79) 

satisfies the triangle inequality and thus introduces a metric. This very quantity can be used 
to bound the trace distance -Dtr(pi,P2) = ^Trjpi — P2I < C{pi,p2) from above [32111]. Note 
that the Bures distance, Bures length and root infidelity are functions of the same quantity, 
so they generate the same topology. 

7.3 Modified Bures length 

In analogy to ([75| and ([76]) one may ask whether 

Dg iPi , P2 ) = v^2-2yG(^^. (80) 

and 

Dc{pi,P2) ■■= arccosv/G(/9i,p2) (81) 

define distances. This is obvious for = 2 because F ^ G. The situation changes for > 3, 
for which the Dp and Dq do differ and only F < G is valid. 

We do not know, whether Dq and D'q are distances. However, it can be proved that a 
direct analogue of the root infidelity (ff^l) 

G'(pi, p2) - V^-G{puP2) (82) 

is a genuine distance. The same is true for the modified Bures length, 

D'aApi^P^) = arccosG(pi,p2)- (83) 

Proof. Let us call C the direct sum of the real linear space of all Hermitian operators 
and the 1-dimensional space of real numbers. Its elements are {H, x}, H Hermitian, x a real 
number. £ becomes Euclidean (i.e. a real Hilbert space) by defining the scalar product 

{{Hi,xi},{H2,X2}) = tiHiH2+xiX2. (84) 

Let us denote by B{£) the unit ball of £ and by S{C) the unit sphere. Our proof rests on 
the embedding of the Hermitian operators 



Bn = {H\ trH = l,trH^ < 1} 



(85) 
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into S{C) by 

H^^h:= {h, Vl-trm} ■ (86) 
Clearly, {£,h,£,h) — 1, and from ([84]) we get 

{Ch,^h')^G{H,H'). (87) 



Now it is obvious that \/2 — 2G is the Euclidean distance between and provided H 
and H' belong to B^. Because the density operators form a subset of Sat, (|82p is a distance. 

From ((87|) we get G(i7, = cos a, where a is the angle from which and are seen 
from the center of the ball B{C). Thus, arccosG = a and, in particular, ([83|) is a distance. 
□ . 

Let us now return to the two conditions of (|85p . They are equivalent with 
En = {H\ iiH = 1, tr(iJ - — I)^ < 



N ' - N 

and they describe the smallest ball containing the state space. Em is an afHne translate 
by 1/N of the generalized Bloch-ball, [8j. Bm is centered Si% A — N~^\ and is of radius 
^{N-l)/N. 

Above we have embedded Em by the map ((86|) into the sphere S{C). Just this gives the 
opportunity to apply Mielnik's definition [40] for a transition probability (he also called it 
afSne ratio) of two extremal states of a compact convex set. In our case the compact convex 
set is B{C) and its extremal part is S[C). At the case at hand, Mielnik's procedure starts 
with first choosing an extremal point ^ e S{C) and selecting all affine functions I satisfying 
— 1 and < Z < 1 on B{£). Any such function can be written 

m = (89) 

a + 2 

with a > and rj ^ C arbitrarily. Now we have to vary over all these affine functions, 

Pm{v,0 min/(77) = min ° ^ , (90) 

( a 2 + a 



to get Mielnik's transition probability 



PMiLv)^ ^ ■ (91) 

Returning to H, H' G B^, we can write 

Pm{H,H') ■.= PM{iH,iH') = ^ + ^(^^'^') (92) 

and, in becoming even more special by choosing two density operators for H and H' in the 
equation above, we arrive at 

a 



Dm{pi,P2) = 2a/1 -pm{pi,P2) = 2 sin-, (93) 



(using 2 cos^ = 1 + cos) and also at 



D'f^i (pi , P2 ) = 2 arccos \Jpm{pi,P2)- (94) 
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8 Concluding remarks 

In this paper we analyzed various bounds for quantum fidelity. Two quantities, we propose to 
call sub- and super-fidelity, posses particularly nice properties. On one hand these quantities 
form universal lower and upper bounds for the fidelity. Moreover, with respect to the tensor 
product they display sub- and super-multiplicativity. 

On the other hand, quantities E and G are much easier to calculate than the original 
fidelity F. To compute any of these bounds it is enough to evaluate three traces only. Thus 
one can expect, the quantities introduced in this paper might become useful for various tasks 
of the theory of quantum information processing. Furthermore, under a realistic assumption 
that several copies of both states are available, it is possible to design a scheme to measure 
experimentally sub- and super-fidelity between arbitrary mixed states. For instance, the 
measurement of super-fidelity is possible if one has three copies of each state. In this paper we 
have worked out concrete schemes of such experiments concerning the super-fidelity between 
any two mixed states representing the polarization of photons. 
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Note added 

After this paper was submitted we learned about a related work by Mendonca et al. [43] 
in which the super-fidelity was independently introduced and was called an 'alternative fi- 
delity' measure. In this valuable work the authors provide an alternative proof of super- 
multiplicativity of G, discuss its relation to the trace distance and analyze the distance G 
induces into the space of mixed quantum states, and prove that G is jointly concave in its 
two arguments. 
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Appendix A Algebraic facts 

In this appendix we collect useful algebraic facts, which are used in the main body of the 
paper. 

Fact 1 (From corollary IX. 5. 3 in [7]) For any positive matrices A and B and every uni- 
tarily invariant norm 1 1 1 - 1 1 1 we have 



where v e [0, 1]. 

Fact 2 (From corollary IX. 5. 4 in ^) For any positive matrices A and B and every uni- 
tarily invariant norm 1 1 1 - 1 1 1 we have 



where v e [0, 1] . 

Next two facts can be found in [16] . 
Fact 3 Matrix AB is similar to matrices V^B\^ and ^/BA^/B. 

Fact 4 For positive matrices A and B matrix AB has positive eigenvalues. 

Fact 5 If pi + P2 + ■ ■ ■ + Pn = 1 and pt > then 




(A.l) 




(A.2) 




(A.3) 



Proposition 4 Let g be defined as 




(A.4) 



For x,y € R" such that 



k 



k 



n 



]^yi, for fc = 1, . . . ,n. 



(A.5) 



i=l 



with equality for k = n, we have 



g{x) < g{y). 



(A.6) 



Proof. We introduce notation 




(A.7) 



Direct computation shows that function g satisfies 



Uigi{u) > W2.92(w) > ■ ■ • > w„.9„(u), 



(A.8) 
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for u e R" such that ui > U2 > ■ ■ ■ > Un > 0. We denote ai = log{xi) and Pi = log(t/i). Note 
that (jA.5|l can be rewritten as 

k k 

'^a,<^P,fork = l,...,n, (A.9) 

1=1 i=l 

with equahty for k = n. We define new function 

h{v)^g{e^\e^\...,e'^-). (A.IO) 
For a given vector u such that Wi > W2 > ■ • • > w„ > 0, and = log(wi) we have 

vi>V2>--->Vn. (A.ll) 

Using (|A.8[) we can write 

e''igi(e"S e''^ . . . , e"") > • • • > e^^ gn{e-\ e^^ . . . , e^"). (A.12) 

Now from above and (jA.lOp we have 

hi{v)>h2{v)>--->hn{v). (A.13) 

Note now that function h satisfies condition from |15[ Theorem 3.6] and thus it is Schur- 
convex, so 

h{a) < h{P). (A.14) 

Using (jA.10[) we can write 

gix) < 9{y)- (A.15) 

Thus the proof is complete. □. 

Fact 6 (Holder's inequality [19j) For a > l,b = a/{a — 1) and positive semidefinite A 
and B we have 

U{AB) < (trA'')i/''(trB'')i/''. (A.16) 
Fact 7 For density matrices A and B we have 



1 - VtrAVtrB^ >\Jl- tr^Vl - trS^. (A.17) 
Proof. This incquafity can be rewritten in equivalent form 



1 - 2VtrA2trB2 + trA^trB"^ > 1 - iiA^ - trB^ + trA^trS^ (A. 18) 

which is equivalent to 

trA^ + trS^ 



Vtryl2trB2 < 1 . (A.19) 

This completes the proof since for any positive numbers the arithmetic mean is always greater 
than or equal to the geometric mean. □. 
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Fact 8 (Maclaurin inequality \42\ p. 5]) For a given matrix A of rank r and with r pos- 
itive eigenvalues we have 

.£p> .Jf^ (A.20) 



\kl V \k+l/ 



for 1 < k < r. 



Appendix B Proof of the lower bound ()10p 

To prove that sub-fidelity E is not larger than fidelity F, let us take a look at equations 
and (j27p in which both quantities are expressed in terms of the second symmetric function. 
We can rewrite the function S2, which forms fidelity, 



S2{VAy^BA^) = A.(v/A1/2b^i/2)Aj-(v/^i/2s^i/2) (B.l) 



J2 ^ k{A^/^BA^/^)J Xj{A^/^BA^/-^) (B.2) 



= Y.;jK{AB)^\,{AB). (B.3) 

The last equality is the consequence of similarity of matrices A'^/'^BA^/'^ and AB. Making 
use of subadditivity of square root we obtain 



S2 



(VAV2bIV2) > fY\,{AB)\j{AB) (B.4) 



= y/s2{AB). (B.5) 

As a consequence we get 



F{A, B) = tvAB + 2s2(V^4V2B^4i72) > trA^ + 2y/7^(AB) = E{A, B). (B.6) 
Appendix C Proof of the lower bound (|lip 

To prove inequality (fTTj) we use Fact [8] (Maclaurin inequality) and obtain 

1/2 , / s . 1/r 



S2 (^^1/2^^1/2^ \ (Vyll/2BAl/2^ 

1 - I TF\ I ■ (C^-l) 



Using Fact [3] we get 

S2 [Va^/^BA^^) > Q (,s,(Vv4i/2iJAi/2) Q [nA.(VAV2i?Ai 



2/r 



/2 



Vi=l 



y (^n (AV2SA1/2) j = a. (ab)^ 



VMAB). 
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Now using (fT5|) we write 



F(A, B) ^ trAB + 2s2 iyA^^^^BA^j > tr AB + r{r - 1) </7jAB). (C.2) 
Appendix D Proofs of Lemmas 

Proof of Lemma [TJ Observe that the matrix p}^^P2Pi^^ is similar to pip2 and thus 



E/, / 1/2 l/2\ > / 1/2 1/2 
Y ^» (Pl P2Pi' j (^p/ p2Pi' 



= X! V^* (^1^2) Aj (P1P2), 

where Ai (A) denotes i^^ eigenvalue of a matrix A. 

Let us define a function g : i?" —f R which acts on a vector x of non-negative numbers 

■■^^y^XtXj. (D.l) 



It allows one to rewrite 



2s2 (^\/pF^P2pP^^ =g(A(pip2)), 



(D.2) 



where X{A) denotes the vector of eigenvalues of A. 
From [18, Theorem 3.3.2 and 3.3.4] we obtain 



Y[kipiP2) <Y[^^{pl)X^ip2) for = 1 , . . . , 71, (D.3) 

1=1 2=1 

with equality for k = n. Making use of Proposition [5] from Appendix A with Xi ~ \i{piP2) 
and Hi = \i{pi)\i{p2) we obtain 

5(A(piP2)) <.g(A(pi)oA(p2)), (D.4) 
where o denotes Hadamard product, 18, Definition 7.5.1]. Now making use of (|D.2p we get 

S2 < S2 (^v^diag(A(pi))diag(A(p2))) • (D.5) 

And thus the proof is complete. □. 

Proof of Lemma [2l For given density matrices pi , p2 with eigenvalues pi , . . . , p„ and 
qi, . . . ,q„ respectively. We denote diagonal matrices with entries on diagonal pi , . . . , p„ and 
qi, . . . ,qn as diag(p), diag(q) respectively. 

Rewriting the second elementary function §2 we obtain 

2s2 (Vdiag(p)diag(g)j = ^ VPj ■ (D-6) 
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On the other hand 

Let us define vectors x,y £ R"^ 

a;ij = y/PiPj{^ - Sij), = y/qiOjil - (D.8) 
where Xij = Xn{i-i)+j- Using Cauchy-Schwarz inequality 

\{x\y)\<^/W)^/My), (D-9) 

we get 

This completes the proof. □. 



Appendix E The case iV = 3 

In this section we are going to study the fidelity of two states pi and p2 in the case where 
the rank r of their product piP2 is not greater than 3. As in Section [2] we will denote 

eigenvalues of ^ P\'^ P2P\'^ by Ai, so eigenvalues of p^"^ p^p^^"^ are given by {Af} and by 
similarity we have that the eigenvalues of p\P2 are also given by {A^}. Since r < 3, not 
more than three eigenvalues of p\p2 are positive, so the third symmetric function ([6]) reads 
s-i{,P\P2) = (AiA2A3)^. This is so for any two states of a qutrit, so for A = 3 one has 
s-i[,P\P2) = det(pip2)- 

Consider now the expression for fidelity (j4]) which can be rewritten with the use of the 
second symmetric function, 

F(pi,p2)=trpip2 + 2s2(\/^FW^)- (E.l) 

The square of the symmetric function presented in the above equation, can be written as 
(X]i<j ^i^i) = Si<j + ^- The reminder R, defined imphcitly by this equation, is equal 
to zero if r < 2 and the sum consists of a single term only. It is difficult to handle R generally. 
But if r = 3 one has 

i?= AiA2(A2A3 + A3Ai) + ... = 2(AiA2A3)(Ai + A2 + As). (E.2) 

However, in the particular case r < 3 discussed here used to ([3]) one has Ai + A2 + A3 = ^J~F 
while A1A2A3 — \fs^yPiP2)- Combining this with (|4]) we get the equation for fidelity satisfied 
for r < 3 

F = trpip2 + S2(piP2) + 2VfVs3(pi^. (E.3) 

In the case r < 2 the third function S3 vanishes, so this equation leads to an expression, 
F = trpi/72 + 2^/ S2{piP2) = tr/3i/32 + 2A1A2, already discussed in Section[2l 
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Another relation for fidelity is due to the fact that an assumption r < 3 implies that 

^ A,A;t = (AiA2A3)(Ar' + A2-1 + Ag-^). (E.4) 

j<k 

Therefore in this case one has 

S2(V^F^^BAV2) = ^/det{AB)F{l/A,l/B). (E.5) 

Lifting for a moment the assumption that the arguments of fidelity have to be normalized, 
we arrive therefore at another equation for fidelity satisfied for r < 3, 

F{A, B) = trAB + 2^ydet{AB)F{l/A, 1/B). (E.6) 



